Search for: All records

Mature faults with large cumulative slip often separate rocks with dissimilar elastic properties and show asymmetric damage distribution. Elastic contrast across such bimaterial faults can significantly modify various aspects of earthquake rupture dynamics, including normal stress variations, rupture propagation direction, distribution of ground motions, and evolution of off‐fault damage. Thus, analyzing elastic contrasts of bimaterial faults is important for understanding earthquake physics and related hazard potential. The effect of elastic contrast between isotropic materials on rupture dynamics is relatively well studied. However, most fault rocks are elastically anisotropic, and little is known about how the anisotropy affects rupture dynamics. We examine microstructures of the Sandhill Corner shear zone, which separates quartzofeldspathic rock and micaceous schist with wider and narrower damage zones, respectively. This shear zone is part of the Norumbega fault system, a Paleozoic, large‐displacement, seismogenic, strike‐slip fault system exhumed from middle crustal depths. We calculate elastic properties and seismic wave speeds of elastically anisotropic rocks from each unit having different proportions of mica grains aligned sub‐parallel to the fault. Our findings show that the horizontally polarized shear wave propagating parallel to the bimaterial fault (with fault‐normal particle motion) is the slowest owing to the fault‐normal compliance and therefore may be important in determining the elastic contrast that affects rupture dynamics in anisotropic media. Following results from subshear rupture propagation models in isotropic media, our results are consistent with ruptures preferentially propagated in the slip direction of the schist, which has the slower horizontal shear wave and larger fault‐normal compliance.

Seismological fracture or breakdown energy represents energy expended in a volume surrounding the advancing rupture front and the slipping fault surface. Estimates are commonly obtained by inverting ground motions and using the results to model slip on the fault surface. However, this practice cannot identify contributions from different energy‐consumption processes, so our understanding of the importance of these processes comes largely from field‐ and laboratory‐based studies. Here, we use garnet fragment size data to estimate surface‐area energy density with distance from the fault core in the damage zone of a deeply exhumed strike‐slip fault/shear zone. Estimated energy densities per fragmentation event range from 2.87 × 10³to 2.72 × 10⁵ J/m³in the outer and inner portions of the dynamic damage zone, respectively, with the dynamic zone being inferred from the fractal dimensions of fragment size distributions and other indicators. Integrating over the ∼105 m width of the dynamic damage zone gives fracture surface‐area energy per unit fault area ranging from a lower bound of 6.63 × 10⁵ J/m²to an upper bound of 1.63 × 10⁷ J/m²per event. This range overlaps with most geological, theoretical, and kinematic slip‐model estimates of energy expenditure in the source volume for earthquakes characterized by seismic moments >10¹⁷ N·m. We employ physics‐based fragmentation models to estimate equivalent tensile strain rates associated with garnet fragmentation, which range from 5.42 × 10²to 1.04 × 10⁴ s⁻¹per earthquake in the outer and inner portions of the dynamic damage zone, respectively. Our results suggest that surface‐energy generation is a nonnegligible component of the earthquake energy budget.